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Abstract 

Using first-principles calculations, we systematically study the dissociation of O2 molecules on 
different ultrathin Pb(lll) films. Based on our previous work revealing the molecular adsorption 
precursor states for O2, we further explore that why there are two nearly degenerate adsorption 
states on Pb(lll) ultrathin films, but no precursor adsorption states exist at all on the Mg(0001) 
and Al(lll) surfaces. And the reason is concluded to be the different surface electronic structures. 
For the O2 dissociation, we consider both the reaction channels from gas-like and molecularly 
adsorbed O2 molecules. We find that the energy barrier for O2 dissociation from the molecular 
adsorption precursor states is always smaller than from O2 gases. The most energetically favorable 
dissociation process is found to be the same on different Pb(lll) films, and the energy barriers are 
found to be modulated by the quantum size effects of Pb(lll) films. 
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I. INTRODUCTION 



The adsorption and dissociation of O2 molecules on metal surfaces are of great impor- 
tance to the subsequent oxidation reactions and to the formation of metal oxides 1 . This is 
especially true for the formation of thin oxide films that have been widely used as catalysts, 
sensors, dielectrics, and corrosion inhibitors 2 . Thus vast studies have been carried out on 
the O2 adsorption and dissociation on metal surfaces. During these studies, the theoret- 
ical ab initio modeling based on the adiabatic approximation has been proved successful 
over a wide range for studying the adsorption and dissociation of O2 on transition metal 
surfaces. By calculating the adiabatic potential energy surface (PES), it has been found 
that O2 molecules will spontaneously dissociate while adsorbing at reactive transition metal 
surfaces like iron (Fe) 3 . For noble transition metals like gold (Au)^, silver (Ag)-, Copper 
(Cu)£, platinum (Pt)^£ and Nickel (Ni)^, the adsorption of 2 turns out to depend on the 
ambient temperature, and both atomic and molecular adsorptions have been observed. Re- 
markably, in all above transition metal systems, the concept of adiabatic calculation works 
very well in explaining and predicting a large amount of physical/chemical phenomena dur- 
ing dissociation process of O2. When the attention is focused on the nontransition metals 
with only sp valence electrons, an uncomfortable gap opens between ab initio prediction 
and experimental observation. The most notable is the long-term enigma of low initial 
sticking probability of thermal O2 molecules at Al(lll), which has been measured by many 
independent experiments 8,9 but cannot be reproduced by adiabatic state-of-the-art density 
functional theory (DFT) calculations^ - —. The central problem is that the adiabatic DFT 
calculations were unable to find any sizeable barriers on the adiabatic PES, which has led to 
speculations that nonadiabatic effects may play an important role in the oxygen dissociation 
process at the Al(lll) surface^ 1 ^"— . Recently, we present a comparative study on the elec- 
tronic structure of an 2 molecule in close to the Be(0001), Mg(0001) and Al(lll) surfaces, 
and find that the triplet state of O2 is influenced by the Mg(0001) and Al(lll) surfaces, but 
not by the Be(0001) surface^. Thus we prove that the adiabatic DFT calculations are still 
reliable to study the O 2 /Be(0001) system, and find out sizable dissociation energy barriers 
for O2 molecules on the Be(0001) surface 1 -. However, in contrast to the systematical re- 
sults on the O2 dissociation on transition metal surfaces, there are still no criteria to judge 
whether an energy barrier is needed or whether a precursor molecular-adsorption state exists 
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for the O2 dissociation on sp metal surfaces. Motivated by these backgrounds, in this paper 
we investigate the dissociation of O2 molecules on another sp metal surface, lead (Pb). We 
find that the adiabatic DFT calculation also leads to reliable results for the 02/Pb(lll) 
system. 

Another motive for us to study the O2 dissociation on the Pb(lll) surface is that atom- 
ically uniform Pb films with thickness from several to tens of monolayers (MLs) have been 
successfully fabricated on both metal and semiconductor substrates 2 ^ 2 ^, which can serve as 
an ideal prototype for experimental researchers to investigate the surface oxidation of Pb. 
However, it is believed that for very thin metal films, the adsorption and dissociation of O2 
molecules becomes more complex because of the two boundaries. According to the "particle- 
in-a-box" model, electrons confined in a uniform film are quantized into discrete energy levels 
along the film-normal direction, forming the so-called quantum well states. Because the elec- 
tron Fermi wavelength of Pb is nearly 4 times the lattice spacing along the [111] direction^ 2 -, 
many properties like thermal stability, surface energy, and work function of the Pb films 
oscillate with thickness with a quasibilayer period^— . Recently, a unique two-step method 
has been brought out to efficiently oxidize the the Pb(lll) surface 2 ^ 1 ^. During this special 
oxidation process, O2 molecules are firstly introduced at a low temperature and adsorb in a 
precursor state, then the Pb(lll) surface is oxidized by a subsequently annealing to room 
temperature 2 ^. Most interestingly, both the oxygen coverage at low temperature and oxide 
coverage after annealing are found to be modulated by the quantum size effect (QSE) of 
ultrathin Pb(lll) films 26 . These studies together with other investigations 28 - clearly prove 
the modulation of surface oxidation by QSE. To understand the complicated experimen- 
tal observations, the systematical investigations on the adsorption and dissociation of O2 
molecules on ultrathin Pb(lll) films are needed. 

In our previous work, we have already studied the adsorption properties for both molec- 
ular O2 and atomic O on different Pb(lll) films 2 ^ 1 ^. It is found that O2 molecules will 
not spontaneously dissociate on Pb(lll) films, and a precursor molecular adsorption state 
exists 2 ^. Here in this work, we further analyze the electronic interactions for the adsorption 
of O2, investigate the dissociation process for O2, and study the QSE on the dissociation 
energy barriers. The rest of the paper organized as follows. In Sec. II, we give details of 
the first-principles total energy calculations, which is followed in Sec. Ill by our analysis for 
the adsorption state of O2 on the Pb(lll) surface. In Sec. IV, we study the O2 dissociation 
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on Pb(lll) films, and discuss the QSE on the corresponding energy barriers. And finally in 
Sec. VI, we give our conclusions. 

II. COMPUTATIONAL METHODS 

The DFT total energy calculations are carried out using the Vienna ab initio simu- 
lation package 31 . The generalized gradient approximation (GGA) of Perdew et al— and 
the projector-augmented-wave (PAW) potentials^ are employed to describe the exchange- 
correlation energy and the electron-ion interaction, respectively. The so-called "repeated 
slab" geometries^ are employed to model the clean Pb(lll) films ranging from 4 to 8 mono- 
layers (ML). A vacuum layer of 20 A is adopted to separate the Pb(lll) slabs in two 
adjacent cells, which is found to be sufficiently convergent from our test calculations. The 
O2 was placed on one side of the slab, namely on the top surface, whereas the bottom layer 
was fixed. All other Pb layers as well as the oxygen atoms are free to relax until the forces 
on them are less than 0.02 eV/A. The plane-wave energy cutoff is set to 400 eV. A Fermi 
broadening 35 of 0.1 eV is chosen to smear the occupation of the bands around Ep by a 
finite- T Fermi function and extrapolating to T=0 K. Integration over the Brillouin zone is 
done using the Monkhorst-Pack scheme^ with 13x13x1 grid points. The calculated lattice 
constant of bulk Pb is 5.03 A, in good agreement with the experimental value of 4.95 A—. 
The binding energy of spin-polarized O2 is calculated to be D=5.78 eV per atom and the 
0-0 bond length is about 1.235 A. These results are typical for well-converged DFT-GGA 
calculations. Compared to the experimental 38 values of 5.12 eV and 1.21 A for O binding 
energy and bonding length, the usual DFT-GGA result always introduces an overestimation, 
which reflects the theoretical deficiency for describing the local orbitals of the oxygen. We 
will consider this overbinding of O2 when drawing any conclusion that may be affected by 
its explicit value. 

III. THE ADSORPTION PROPERTIES OF 2 ON THE PB(lll) SURFACE 

The geometries of ultrathin Pb(lll) films are firstly relaxed before our study for the 
adsorption and dissociation of O2 molecules. It is found that the surface relaxations display 
well-defined QSE™. Generally speaking, atoms at the outmost layer tend to relax inward. 
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In contrast, atoms at the second layer tend to relax outward. And both the contraction for 
the first interlayer and the expansion for the second interlayer exhibit oscillations of bilayer 
period, reflecting the effect of the quantum well states^. 

We then check the reliability of our adiabatic calculations. The main reason for the 
wrong conclusions lead by adiabatic first-principles calculations on the adsorption of O2 on 
sp-metals lies in that the lowest unoccupied molecular orbital (LUMO) state of O2 is always 
aligned with the Fermi level at any distance between the molecule and metal surfaces such 
as Al(lll), allowing a partial filling of the empty molecular orbital 3 -. Due to this unphysi- 
cal (metal-molecule or intramolecule) charge transfer, the adsorption of O2 on the Al(lll) 
surface is calculated to be dissociative without any energy barriers, which contradicts with 
the experimental observations of low initial sticking probability at low temperatures^. In 
our previous work—, we have shown that the unphysical partial filling of the LUMO does 
not happen during the adsorption of 2 on the Be(0001) surface, and thus we obtained 
sizeable energy barriers for the dissociation of O2. Here for the 02/Pb(lll) system, molec- 
ular adsorption precursor states exist, so we can not judge the reliability of our calculations 
by whether a sizable energy barrier exists or not. We calculate and show in Fig. 1(a) the 
electronic density of states (DOS) of an O2 molecule in close to a 5 ML Pb(lll) film. In 
comparison with the DOS of O2 in close to the Mg(0001) and Al(lll) surfaces depicted in 
Figs. 1(b) and (c), we find that the triplet ground state of the O2 molecule is not influenced 
at all in close to the Pb(lll) film, but totally disturbed by the Mg(0001) and Al(lll) sur- 
faces due to the unphysical charge transfers. Therefore, we are sure that the adiabatic DFT 
calculations will not lead to unphysical conclusions for the 02/Pb(lll) system, as well as 
for the O 2 /Be(0001) system. 

To facilitate our discussions, we give some notations at first. Formally, there are four 
high-symmetry sites on the Pb(lll) surface, respectively the top (T), bridge (B), hep (HH), 
and fee (FH) hollow sites. Our previous study has already shown that the high-symmetry 
sites play crucial roles in the adsorption O2 molecules on Pb(lll) films^. Therefore we do 
not consider other surface sites any more. At each surface site, an O2 molecule has three 
different high-symmetry orientations, respectively along the x (i.e., [Oil]), y (i.e., [211]), and 
z (i.e., [Ill]) directions. Therefore, we will use T-x, y, z, B-x, y, z, HH-x, y, z, and FH-x, y, z 
to respectively represent the 12 high symmetry adsorption channels for O2. 

After geometry optimizations by initially setting the O2 molecule at 4.0 A along each 
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adsorption channel, we have got 6 molecular adsorption states for O^. By using numbers 
1, 2, and 3 to represent the O2 orientation along the x, y, and z directions, we name the 
six adsorption states respectively as HH1, HH2, HH3, FH1, FH2, and FH3. The atomic 
geometries for the HH1, HH2 and HH3 states are listed in Figs. 2(b), (c) and (d) respectively. 
The geometries for the FH1, FH2 and FH3 states are very similar to HH1, HH2 and HH3 
states, with the adsorbed oxygen atoms at fee hollow sites instead of hep hollow sites. Among 
the six precursor states, HH1 and HH2 are the most stable and two nearly degenerate states. 
In our previous work, we have already discussed in detail the interactions between oxygen 
and Pb for the HH1 stated Here we further analyze the electronic interactions for the HH2 
state to discuss why the two adsorption states with different O2 orientations (i.e. HH1 and 
HH2 states) are almost degenerate. 

Figure 3 shows the calculated difference charge density Ap(r) for the HH2 state on the 5 
ML Pb(lll) film, which is obtained by subtracting the electron densities of noninteracting 
component systems, p Pb ( m )(r) + p° 2 (r), from the density p(r) of the HH2 state, while 
retaining the atomic positions of the component systems at the same location as in HH2. 
A positive Ap(r) here represents charge accumulation, whereas a negative Ap(r) represents 
charge depletion. Figure 3(a) clearly shows a charge transfer from Pb to O2, which is 
consistent with the electronegativity difference between O and Pb atoms, and same with 
what we found for the HH1 state 29 . The contour plot for Ap(r) in the plane containing the 
two O atoms and parallel to the Pb surface is shown Fig. 3(b), we can see that there is 
an electron depletion area in the middle of the two oxygen atoms, indicating the donation 
of electrons from bonding orbitals of O2 to Pb. We can also see two electron accumulation 
areas around the two oxygen atoms indicating the back donation of electrons from Pb to 
the 7T* antibonding orbital of O2. This interaction mechanism of CVto-surface donation and 
surface-to-02 vr* back donation is just the same with that in the HH1 stated. So it is no 
wonder that the two adsorption states are nearly degenerate and have similar adsorption 
energies. 

In comparison with the Be(0001), Mg(0001) and Al(lll) surfaces, Pb(lll) is unique in 
that a parallel O2 molecularly adsorbs only on it. This is because of the different surface 
electronic structures. As shown in Fig. 4, the p electronic states around the Fermi energy 
strongly hybridize with s states for both the Mg and Al surfaces (also for the Be(0001) 
surface^), but does not hybridize at all for the Pb(lll) surface. As we have seen~, the 
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interaction between O2 and the Pb(lll) surface is mainly contributed by the electronic 
hybridizations between the frontier molecular orbitals (the orbitals nearest to the Fermi 
energy) of O2 and electronic states of Pb. According to the frontier orbital theory, the 
interactions between O2 and other metals have similar mechanisms. However, the frontier 
orbitals of 2 (i.e. a p , ir p , and 7r*) are composed of p electrons of oxygen, thus they are more 
compatible with p electronic states than sp hybridized electronic states. This compatibility 
between frontier orbitals of O2 and p states of Pb finally leads to the molecular adsorption 
precursor states. 

IV. DISSOCIATION OF 2 MOLECULES 

After the studies and discussions for the adsorption on different Pb(lll) films, we now 
investigate the 2 dissociation. Because of the existence of precursor states, there are two 
different kinds of dissociation processes for O2, respectively from O2 gases and from the 
molecular adsorption states. Because the adiabatic calculations have been proved reliable 
for the 02/Pb(lll) system, we will calculate the adiabatic one- dimensional (ID) and two- 
dimensional (2D) potential energy surface (PES) cuts to evaluate the dissociation energy 
barriers for O2 molecules. 

The ID PES cuts for an O2 molecule along the T-x, y and B-x, y channels are show in 
Fig. 5, together with the 0-0 bond length as functions of the O2 height h from the 5 ML 
Pb(lll) film surface. The adsorption from O2 gas along the FH- and B\B\-x,y channels will 
evolve into the molecular adsorption precursor states. Clearly, the dissociative adsorption of 
O2 along the T-x, y and B-y channels are direct and activated, but no dissociative adsorption 
is found for an O2 molecules along the B-x channel. This is different from the O2/Be(0001) 
system, where the O2 dissociation along B-x channel also occurs^. Among the calculated 
channels, the lowest dissociation energy barrier is found to be 0.36 eV along the B-y channel, 
with the 0-0 bond length and O2 height at the transition state to be 1.33 and 2.20A. The 
ID PES cuts along the T-x, y channels have similar distributions, and the corresponding 
energy barriers are 0.68 and 0.59 eV. 

The electronic interactions are then analyzed for the dissociation process from O2 gases. 
Figure 6 shows the PDOS for an O2 molecule at the initial, transition and final states of 
the dissociation path along the B-y channel. We can see that at the initial state, the O2 
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molecule keeps its triplet electronic structure. As we have discussed, the LUMO of O2 is 
above the Fermi energy, proving the reliability of our adiabatic DFT calculations. The 
highest occupied molecular orbital (HOMO) and LUMO are the spin-up and spin-down 
antibonding 7r* molecular orbitals, respectively. Besides of HOMO and LUMO, p electrons 
of O2 make up of two bonding orbitals, a p and n p . For both spins, the a p orbitals are lower 
in energy than the n p orbitals. When the O2 molecule evolves into the transition state, 
electrons transfer from the a p and ir p molecular orbitals to Pb through hybridizations. As 
shown in Fig. 6(b), the a p and 7r p peaks are diminished. We can also see from Fig. 6(b) that 
the HOMO and LUMO of O2 are broadened through hybridizations with electronic states of 
Pb. During the electronic hybridizations, the spin splittings for all the O2 molecular orbitals 
are all decreased. In the final dissociated adsorption state, the total spin of the adsorption 
system becomes zero. We can see from Fig. 6(c) that the molecular orbitals of O2 no longer 
exist after electronic hybridizations. As discussed in our previous work—, the interactions 
between O and Pb atoms display a mixed ionic/covalent character. 

The dissociation of O2 from the molecular adsorption precursor states are more compli- 
cated than from O2 gases because the formation and breaking of O-Pb bonds during the 
dissociation of O2 also need to be considered. Detailed information on the atomic adsorp- 
tion of oxygen on the Pb(lll) surface can help us to simplify the problem. In our previous 
study 30 , we have found that the O atoms choose to adsorb at hollow sites on the Pb(lll) 
surface. Thus we only need to consider the dissociation of an adsorbed O2 molecule into 
different hollow sites. For an O2 molecule in the HH1 or HH2 state, there are two different 
dissociation paths, respectively to the atomic adsorption of oxygen at two neighboring fee 
and hep hollow sites and to the atomic adsorption of oxygen at two fee hollow sites. The 
calculated 2D PES cuts for the adsorbed O2 molecules in HH1 and HH2 states to dissociate 
into two neighboring hollow sites are shown in Figs. 7(a) and (b), while the calculated 2D 
PES cuts for O2 to dissociate into two fee hollow sites are also shown in Figs. 7(c) and 
(d). The calculated energy barriers for these four paths are respectively 0.22, 0.26, 0.59 and 
0.37 eV. These values indicate that the most energetically favorable dissociation path for an 
adsorbed O2 molecule is from the HH1 adsorption state into the atomic adsorption of oxygen 
at two neighboring fee and hep hollow sites. The saddle point in Fig. 7(a) corresponds to 
the adsorption structure of two oxygen atoms near two neighboring bridge sites, with their 
horizontal distance to be 1.40 A. The smallest energy barriers for an O2 molecule in the 
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FH1 and FH2 states to dissociate are both found to be 0.25 eV, which are larger than the 
smallest energy barrier from the HH1 state of 0.22 eV. 

The electronic properties are then analyzed for the dissociation path shown in Fig. 7(a), 
whose energy barrier is the smallest during the dissociation from molecular adsorption pre- 
cursor states. The PDOS around the two oxygen atoms are calculated for the corresponding 
initial, transition and final states. As shown in Fig. 8(a), the it* antibonding orbital of O2 
are broadened and shifted below the Fermi energy in the HH1 adsorption state, through 
hybridizations with p electronic states of Pb. At the transition state of the dissociation 
process, the energy interval between the a s and a* orbitals is decreased from 6.58 to 6.24 
eV, and the oxygen 2p states hybridize in a different way with electronic states of Pb from 
that in the HH1 adsorption state. At the end of the dissociation process, the two oxygen 
atoms relax into two neighboring hollow sites, and the molecular orbitals of O2 all disappear 
after electronic hybridizations. 

After the systematic studies on the O2 dissociation in different ways, we reveal that the 
most energetically favorable one is from the HH1 adsorption state into two neighboring 
hollow sites. It is always believed that because of the confinement of the two boundaries 
in the vertical direction, the quantum behavior of electrons will cause a number of unique 
properties for thin metal films that could be superior to their bulk counterpart. In particular, 
quantum size effects on surface reactivity have been studied extensively for applications in 
catalysis, corrosion, and gas sensing. Here we further consider the dissociation path for 
O2 molecules on different Pb(lll) films, including both the dissociation along the T- and 
B-a;, y channels from O2 gases, and the dissociation from molecular adsorption states. After 
systematic calculations, we find that the most energetically favorable dissociation paths 
are the same on different Pb(lll) films, i.e. from the HH1 adsorption state to the atomic 
adsorption of two oxygen atoms at two neighboring hollow sites. Figure 9 shows the 2D 
PES cuts for the O2 dissociation on 4~7 monolayers Pb(lll) films. The calculated energy 
barriers are respectively 0.26, 0.22, 0.26 and 0.24 eV. Since many properties of Pb(lll) 
films such as work function and surface energy show bilayer oscillation behaviors because of 
QSE 25 . The even-odd oscillation of the dissociation energy barrier proves the modulation of 
QSE on the O2 dissociation on Pb(lll) films. 
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V. CONCLUSION 



In summary, we have investigated the adsorption and dissociation of O2 molecules on 
different Pb(lll) films. The electronic states of an 2 molecule in close to the Pb films are 
firstly analyzed and imply that the adiabatic calculations are reliable. For the molecular 
adsorption precursor states, we find that the electronic interactions are very similar in the 
HH1 and HH2 states, and they are thus two nearly degenerate states. Besides, we further 
point out that no s-p hybridizations in Pb(lll) films is the main reason for the existence 
of molecular adsorption precursor states. Based on our adiabatic PES calculations, we find 
that the dissociation of 2 from the molecular adsorption precursor states have lower energy 
barriers than from 2 gases. The most energetically favorable dissociation paths on different 
Pb(lll) films are found to be the same, which is from the HH1 adsorption state to the atomic 
adsorption of two oxygen atoms at two neighboring hollow sites. The energy barriers of 2 
on different Pb(lll) films are found to be modulated by the QSE of Pb(lll) films, and 
show a bilayer oscillation behavior. Our study further complete the theoretical picture for 
the adsorption and dissociation of 2 molecules on different Pb(lll) films. 
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List of captions 

Fig.l (color online) The electronic density of states for an 2 molecule in close to the 
5 ML Pb(lll) film (a), the Mg(0001) surface (b), and the Al(lll) surface (c). The Fermi 
energies are all set to 0. 

Fig. 2 (color online) (a) The considered twelve high-symmetry adsorption channels 
for 2 on the Pb(lll) surface, (b), (c) and (d) The three molecular adsorption states 
for 2 at the surface hep hollow site. Red and grey balls represent O and Pb atoms. 
There is a Pb atom in the second layer below the HH site to make it distinct from the FH site. 

Fig. 3 (color online) (a) Difference charge density for the HH2 adsorption state of 2 on 
the 5 ML Pb film. Red and gray atoms represent O and surface Pb. Regions of electrons 
accumulation/depletion are displayed in blue/yellow, respectively and the isosurface value 
is ±0.02 e/A. (b) Contour plots of the difference charge density in the plane that contains 
the two oxygen atoms and is parallel to the Pb(lll) surface. 

Fig. 4 (color online) The projected density of states for the 5 ML Pb(lll) film (a), 
Mg(0001) (b) and Al(lll) (c) surfaces. The Fermi energies are all set to 0. 

Fig. 5 (color online) One-dimensional cuts of the potential energy surfaces and the 
corresponding 0-0 bond lengths (d -o) as functions of the 2 distance h from the 5 ML 
Pb(lll) film for eight different dissociative channels. The inset in each panel indicates the 
initial (red circles) and final (green circles) atomic positions of the two oxygen atoms. 

Fig. 6 (color online) The projected density of states around the two oxygen atoms at the 
initial (a), transition (b) and final states (c) along the B-y adsorption channel. The Fermi 
energies are set to 0. 

Fig. 7 (color online) Color-filled contour plots of the potential energy surface for 2 
dissociation on the 5 ML Pb(lll) film from HH1 (a) and HH2 (b) states to the atomic 
adsorption of two oxygen atoms at neighboring hep and fee hollow sites, and from HH1 (c) 
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and HH2 (d) states to the atomic adsorption at two fee hollow sites, as functions of the 2 
bond length do~o and distance h (from the surface). 

Fig. 8 (color online) The projected density of states around the two oxygen atoms at the 
initial (a), transition (b) and final states (c) of the dissociation path from the HH1 state to 
the atomic adsorption of two oxygen atoms at neighboring hep and fee hollow sites. The 
Fermi energies are set to 0. 

Fig. 9 (color online) Color-filled contour plots of the potential energy surface for 2 
dissociation on the 4 (a), 5 (b), 6 (c) and 7 ML (d) Pb(lll) films from the HH1 state 
to the atomic adsorption of two oxygen atoms at neighboring hep and fee hollow sites, as 
functions of the 2 bond length cfo-o and distance h (from the surface). 
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